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I. INTRODUCTION 


Early radar and communication systems used antenna arrays 
consisting of a combination of separate radiating elements. 
These systems date back to the turn of the century [Ref. 1: p. 
11-2]. The antenna radiation characteristics are determined by 
the position of the radiating elements and the amplitude and 
phase excitation relationships. As the operating frequency of 
radar systems increased (decreased wavelength), array antenna 
radars where replaced by simpler parabolic antennas. With the 
advent of printed circuit technology, electronically 
controlled phase shifters have become more compact and low 
cost. Consequently attention has been directed to array 
antennas again. The aperture excitation may now be modulated 
by controlling the individual elements to give beams that are 
scanned electronically [Ref. 1: p. 11-2]. An electronically 
steered array antenna can track a large number of targets, 
illuminate targets with RF energy to guide missiles toward 
targets, perform a complete hemispherical search with 
automatic target selection, and perform multiple functions in 
communication systems [Ref. 1: p. 11-2]. 

The performance of a phased array antenna depends upon the 
phase and amplitude distribution across the aperture. A 


continuous’) linear phase distribution is required for an 


efficient focussed beam, but practical phase shifters are 
digital devices and can only provide an approximation to a 
linear phase distribution. Consequently some type of roundoff 
criteria must be established. However the method of roundoff 
affects the radiation pattern characteristics such as beam 
location and sidelobe level. Accurate target tracking requires 
that a radar have a small beam pointing error. Low sidelobes 
are also desirable to prevent jamming and the illumination of 
clutter. Therefore the goal is to select a roundoff criterion 
that provides a phase distribution across the aperture to 
minimize both the beam shift and the sidelobe levels, while 
Simultaneously maximizing the gain. 

This paper examines the performance degradation of a linear 
array due to the quantization of phase, which is present in 
all practical electronically scanned arrays. Analysis was 
restricted to linear arrays so that the simulations were 
computationally fast. The same techniques can be extended to 
two and three dimensional arrays. The performance of a 
monopulse antenna were simulated on the computer for a linear 
array using the Bayliss distribution. A Bayliss distribution 
is a low sidelobe difference beam distribution commonly used 
in radar. An ideal difference beam has a perfect null at the 
scan angle for a continuous linear phase distribution. When 
phase and amplitude errors are present, a perfect null is not 
possible. In practice a good null is about 35 dB below the sum 


beam peak {Ret. 0p.) 2.272 


The size of the arrays studied was be limited to 26, 50, 
and 76 radiating elements using 3, 4, and 5 bit phase 
shifters. A phase bit 1s a measure of how accurately a phase 
shifter can represent a linear phase by a step approximation; 
the more bits, the closer the phase steps. Although a large 
number of bits reproduce a linear phase better, they add cost 
and complexity to the antenna. Thus there is a tradeoff 
between the bitsize and antenna performance. In addition to 
the phase bitsize, the method of roundoff also affects the 
antenna performance. Four different methods of phase roundoff 
were examined to access their impact on antenna performance. 
The methods examined are referred to as regular roundoff, 
weighted random roundoff, running sum roundoff and symmetric 
running sum roundoff. The first two are in common use, but the 
third and fourth are new methods examined in this paper. It 
was demonstrated that the latter two will have the minimum 
beam shift of the four roundoff methods, without significantly 
degrading the other pattern properties. The antenna system 
parameters evaluated include the beam pointing error, maximum 
Sidelobe level, and the null depth of the difference beam. 

This paper is divided into the following sections: basic 
antenna array theory, phase shifters, phase quantization 


algorithms, and analysis of the data. 


II. PHASE SHIFTERS: AN OVERVIEW 

The evolution of sophisticated phased array antennas is due 
to two factors: the development of electronically controlled 
phase shifters and the integration of computers into the 
system. Electronically controlled phase shifters allow the 
beam to be redirected in less than a few microseconds and the 
computer can process a large amount of data. The principle 
advantage of these two factors is the minimum amount of time 
spent moving a beam between targets and making decisions about 
detection and tracking [Ref. 1: p. 11-2]. The performance of 
a phased array antenna is highly dependent upon its phase 
shifters. An N bit phase shifter will have a most significant 
bit of nm radians (180 degrees) and a least significant bit of 
2m/2% radians. For best performance, a large number of bits 
should be used in the phase shifters, while for minimum cost 
and complexity, a small number of bits is desired. Since most 
feed methods require a phase shifter at every element of the 
array, the performance, size and cost of the antenna system 
are strongly affected by the phase shifter {Ref. 1: 12-3,4]. 

Several parameters of phase shifter design impact the 
performance of the phased array radar. The issues of concern 
are insertion loss, switching time, drive power, phase error, 
and physical size. The insertion loss should be as low as 


possible to reduce the loss of power during transmission and 


lower the signal-to-noise ratio at the receiver. A low 
insertion loss also reduces phase shifter heating due to ohmic 
losses. If amplifiers are used between the phase shifter and 
the radiating element, the power capacity and insertion loss 
become insignificant. The phase shifter switching time should 
be as short as possible. A long switching time increases the 
minimum radar range when non-reciprocal phase shifters are 
used and when a burst of pulses are transmitted in different 
@erections (Ref. 1: p. 12-3,4]. 

The drive power should be as small as possible. A large 
amount of drive power generates heat and may require power 
Supplies that are too large for the given system. They also 
require expensive driver circuit components [Ref. 1: p. 12-3]. 
The phase error should also be as small as possible. It should 
not reduce the antenna gain substantially on transmitting or 
raise sidelobes on receiving. One cause of phase error is the 
size of the least significant bit of digital phase shifters. 
Others are due to the tolerances in the phase shifters and 
driver. The maximum acceptable physical size is limited by the 
element separation (usually A/2 by A/2) in order to be 
packaged behind each radiating element [{Ref. 1: p. 12-3,4]. 

Today there are two types of phase shifters used in radar 
systems: ferrite and diode. Ferrite phase shifters are analog 
devices, although in practice they are digitally controlled. 
The first ferrite phase shifter employed in an operational 


phased array was the Reggia-Spencer phase shifter. The 


development of a latching toroidal ferrite phase shifter, 
which employs a transverse magnetizing field, made significant 
improvements. The driving power was reduced and the driver 
circuitry was simplified. Also, the switching time was reduced 
from about a millisecond to a few microseconds. Ferrite phase 
shifters are impractical in radar systems operating at L-band 
(1 to 2 Ghz) frequencies and below since the magnetizing field 
1S operating below resonance. This is due to available 
materials, which have a lower limit of saturation moment of 
approximately 200 gauss. At L band, a large saturation moment 
produces peak power limiting at a low power level, in a 
reasonably compact phase shifter. In addition, a low moment 
material is very temperature sensitive [Ref. 1: p. 12-4,5]. 
Diode phase shifters were developed with the recognition 
of the importance of the pn diode in radar systems. As 
fabrication techniques improved the diode breakdown voltage, 
cutoff frequencies, and reliability increased. Diode phase 
shifters are temperature insensitive over a wide range and do 
not require air or water cooling, as do the ferrite phase 
shifters. Use of diode phase shifters in radar systems is 
practical at all operating frequencies. At lower frequencies 
(below L band) diode phase shifters perform better since the 
semiconductor loss decreases as frequency decreases. This 
decrease can be traded for high peak power performance by 
adjusting the impedance of the transmission line on which the 


diodes are mounted [Ref. 1: p. 12-4,5]. 


At frequencies at or above L band, the selection of the 
phase shifter type is more difficult, because of high diode 
phase shifter losses and the difficulty in maintaining 
manufacturing tolerances. Table 1 summarizes the performance 


of the phase shifters at these frequencies [Ref. 1: p. 12-5]. 


TABLE 1: SUMMARY OF PERFORMANCE OF TYPICAL FERRITE AND DIODE 
PHASE SHIFTERS 


INSERTION LOSS : Lt 2202 dB;, 
increases with 
the number of 


bits 


SWITCHING TIME 2.055.0 sec 50 nsec - 2 — 
DRIVE POWER 2 W AT 1,000 pps gO ots 5 W 


TRANSMITTED POWER 75 KW PEAK, 10 kW PEAK, 
400 W AVERAGE 200 W AVERAGE 





III. ANTENNA ARRAY THEORY 

In the following discussion an antenna array refers to a 
group of similar elements arranged in an arbitrary geometric 
configuration, and excited with a prescribed phase and 
amplitude relation to give a desired radiation pattern. The 
radiation characteristics most commonly examined are the 
direction and width of the main beam, sidelobe levels, and 
null depth or main lobe level. In this chapter the 
fundamental theory and characteristics of linear arrays will 
be examined. 

The total radiation pattern is the sum of the contributions 
from the individual elements. For a linear array having N 


radiators, with element n excited by amplitude, a., and phase, 


nf 


y, and inter-element spacing 5, the) Padtatton pateemnur. 


gf Gea +¥,)) £(8) (III-1) 
where £(8) is the radiation pattern of a single element [Ref. 
1: p. 11-9]. A typical linear array is shown in Figure 1. The 
term J includes the inter-element phase to scan the beanies 
well as any phase errors that may be present. For convenience, 
the scanning phase of 2ns sin(Q,)/A radians per element will 
be separated out so that J is entirely due to the phase error 


(Ref 1. : p. 11-9], 


2 Ebr Neen . . Noi -¢ on . - 
a s(sin (6) -sin (89) ) j(—— ns (sin (6) -sin (65) +,) 
E(0)=f(0)e° 4 ? Yiae 7 <. 





(III-2) 


n=0 


The leading exponential references the phase to the center of 
the array and will be dropped since it represents a phase 


shift common to all elements and does not affect the antenna 


paeclern. 





Figure 1. Linear Array with N Radiators 
Uniformly Spaced by a Distance s 


If there are no phase errors and all amplitudes are equal 
(Wy =0 and eel), then equation III-2 can be reduced to a 


closed form expression, [Ref 1: p. 11-9} 


sin [Nx (sin (8) +sin (8) ) ] 


(III-3) 
Nsin{nx 


(sin (8) +sin(6,) )] 


mlm l~lu 


Equations III-2 and III-3 are the fundamental response 
equations for a scanned array system. As a typical example, 


Figure 2 shows the radiation pattern calculated from equation 
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Figure 2. A Continuous Linear Phase Taylor Distribution 





III-2 , with 50 elements and s = -3A. In this casevally Wome 
zero and the a are chosen te give a 925 dB icy ice 


distribution. The radiation pattern will have only one major 


10 


lobe and grating lobe maxima will not occur for values of -n/2 
< 6 < n/2 as long as [Ref. 1: p. 11-11] 


+= (sin (8) -sin(6,)) <#n . (III—4) 
The beam of an antenna points in the direction normal to the 
phase front. In phased arrays, the phase front is adjusted by 


controlling the phase of each radiating element, as shown in 


Figure 3. 
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Figure 3. Phased Array 





A. RADIATING ELEMENTS 

The most common type of radiating elements used in array 
antennas are slots, dipoles, horns, and more recently 
conductive patches. Since the size of the radiating element is 
constrained to fit in the array geometry, the area of the 
elements are limited to approximately X/4. The predominate 
electrical factors in choosing the radiating element in an 


array environment is the impedance and the pattern of a 


de 


radiator. Another factor that influences the type of element 
is the feeding method. The radiating element should be chosen 
so that the antenna aperture and feed system can be packaged 
to satisfy the size and weight requirements of the antenna. 
For example, if the feed system contains’ stripline 
components, a stripline dipole would be the optimal choice. On 
the other hand, if waveguide feed devices are used, open-ended 
waveguides or small horn antennas are more convenient [Ref 1: 
p. Jassie 

Most well matched antenna elements have element patterns 
that can be approximated by £(@) = cos(@), at least for angles 
near broadside (normal to the antenna aperture). In this 
analysis only the array factor is considered; the element 
pattern will be neglected (£(8) = 1). This simplification will 
not significantly affect the results because the scan angles 
of interest are limited to 0, = + 60°, and comparisons are 


referenced to an error-free array also with £(0) = 1. 
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IV. PHASE QUANTIZATION ALGORITHMS 
In order to steer a focused beam to a certain direction, 
a linear phase must be set up across the aperture according to 
equation III-2. However, most phase shifters are digital 
devices and the allowable phase values are limited by the bit 
size. Therefore, in most instances the phase at each element 
must be quantized. For an bit phase shifter the bit size is 


determined by 


Bit Size=25° | (IV-1) 
Qn 


Each allowable phase shifter value is referred to as a phase 


state. The phase states are given by 


Ope Sins ze) 1=1,2)onn.; 27) (Iv-2) 


For example, a 4 bit device has a bit size of 360°/2*=22.5° and 
16 phase states. The difference between the desired phase and 
the quantized phase is denoted the quantization error. 

In this paper four different methods of phase quantization 
were used to examine the phased array antenna performance. 
They are regular roundoff, weighted random roundoff, running 
sum roundoff and symmetric running sum roundoff. Each of these 
methods of phase quantization is discussed in detail in the 


following sections. 


V2 


To examine the impact of quantization error on the 


difference beam radiation pattern consider equation III-2: 


N-1 nye . : 
at sn (sin(6) -sin(6,)) +¥,) 
E(6)=Ya,e * / (IV-3) 
n=0 


where a, is the Bayliss coefficient and WY is the quantization 
error at the nth element. Equation IV-3 assumes a line source 
and that the only errors present are those due to the phase 


quantization. Therefore at any scan angle 


N-1i 
E(6,) = ae?” (IV-4) 
n=0 


For a difference beam the coefficients a. will be positives 
one half of the elements (1 < n < N/2 , N even) and negative 
for the other half (N/2 + 1 <n < N). Thus, if a null is to 
occur at 8, then |E(6,)| +0, [Ref. 2: p. 2] and the weighted 


sum of the quantization errors must be Zero. 


A. REGULAR ROUNDOFF 

Regular roundoff is the most straightforward of the four 
methods considered. It involves calculating the desired phase 
at each radiating element and then rounding the phase up or 
down to the closest allowable state of the phase shifter. This 
method of phase estimation results in a periodic quantization 


error, the period of which varies with the bitsize and scan 
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angle. A typical desired versus quantized phase angle for 
regular roundoff is shown in Figure 4. This type of phase 
error causes large sidelobes and potentially large beam 


pointing errors. 
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Figure 4. Desired (0) and Regular Roundoff (*) Phase Angle 
Distribution for a 26 Element Array, 3 Bit Phase Shifter 


B. WEIGHTED RANDOM ROUNDOFF 

This method rounds up or down randomly such that the 
expected value of the phase at any element is the desired 
value. This destroys the periodic quantization error, but also 


raises the average sidelobe level. 
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This method is implemented by choosing a random number 
uniformly distributed between 0 and 1 and assigning it to an 
element of the array. The algorithm combines the random number 
and the desired phase to determine the phase state at the 
Given radiating element as follows. At each element the 
algorithm subtracts the closest phase state from the desired 
phase (which gives the regular roundoff quantization error) 
and divides the difference by the bitsize. The algorithm then 
adds this to the random number chosen for this particular 
element. If the sum is greater than one, the phase at the 
radiating element is rounded up, otherwise the phase is 
rounded down. Using this method, the expected value of the 
quantization error tend towards zero as the number of the 
array elements increases. A typical desired versus weighted 
random roundoff algorithm phase angle, wW, is shown in Figure 


ee 


C. RUNNING SUM ROUNDOFF 

This method chooses the phase state at the nth radiating 
element as the one which keeps the magnitude of the cumulative 
sum of the quantization errors the smallest. This minimizes 
the cumulative quantization error at any point along the 
array. For instance, at element n the algorithm compares the 
cumulative quantization error that results from either 
rounding up or down at the element. The algorithm then chooses 


the phase state that minimizes this total quantization error. 
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This condition will give |E(8,)| > 0 as the number of elements 
is increased. A typical desired versus running sum roundoff 


algorithm phase angle, wW, is shown in Figure 6. 


D. SYMMETRIC RUNNING SUM ROUNDOFF 


This method is related to the running sum roundoff just 
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Figure 5. Desired (0) and Weighted Random Roundoff (*) Phase 
Angle Distribution for a 26 Element Array, 3 Bit Phase Shifter 


described, but takes advantage of the phase symmetry across 


the array. If the phase reference is defined at the center of 
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the array, then elements located at equal distances from the 


center will have equal but opposite phase values. For one half 
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Figure 6. Desired (0) and Running Sum Roundoff Phase Angle 
Distribution for a 26 Element Array, 3 Bit Phase Shifter 


of the array, the running sum roundoff algorithm can be used 
to determine the quantized phase. The phase state for the 
corresponding elements on the second half of the array will be 
set so that the quantization error is equal and opposite. 
The rational for this approach is based upon equation IV- 
4, which for a sum beam with a symmetric amplitude 


distribution becomes, 
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N ; 
E(6,)=>)) a,e’* 
n=1 
N (IV-5) 
2 , | 
=) a,(e) M+e) WI) 
n=] 


For small errors, 


N 


2 
: IV-6 
E(6,)=)> an(2+5 +5 Worn) 
n=1 
iis choosing Wf, = - will force the error in each term of 


the sum to zero. In practice this is not possible because the 
cosine term from the exponential, which leads to the 2, is 
Significantly in error. Therefore the peak of the beam is 
somewhat less than 2 a,. This method is included in the study 
to determine whether there is a symmetric property that can be 
exploited. A typical desired versus symmetric running sum 


roundoff algorithm phase angle, WwW, is shown in Figure 7. 


E. QUANTIZATION ERROR AND NULL DEPTH AND LOCATION 

In the case of a difference beam, if the roundoff 
algorithm could set up a linear phase across the aperture at 
a scan angle, then there would be a perfect null at that angle 
(1.e. Zero pointing error and zero gain). Since the bitsize of 


the phase shifter is limited, applying a roundoff algorithm to 


ro 


an array fills the null and introduces a pointing erromeaea- 


shown in Figure 8. The phase error causes the null to be 
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Figure 7. Desired (o) and Symmetric Running Sum Roundoff Phase 
Angle Distribution for a 26 Element Array, 3 Bit Phase Shifter 


rounded instead of a sharp notch, and moves the minimum value 
away from @,. Both the location and depth of the nuligarm 
important. If a null is sharp but located at the wrong angle, 
a tracking error will be introduced in a radar system. If the 


null is too shallow (rounded), then the radar will not be able 
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to position the minimum on the target, and again a tracking 


error will be introduced. 


PERFECT versus ROUNDOFF ALGORITHM NULL 
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Figure 8. Perfect Bayliss Difference Beam Null versus Roundoff 
Algorithm Null 
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V. ANALYSIS OF THE DATA 


In this chapter, data from the different phased array 
antenna configurations will be examined. Comparisons will 
examine the effects on the null depth, sidelobe levels, and 
pointing error as a function of the number of elements, 
roundoff algorithm, and the number of bits. Only one of the 


these will be changed in any given comparison. 


A. COMPARISON BY ROUNDOFF ALGORITHMS 

In this section, typical data is presented for an array 
of 50 elements using the four roundoff algorithms with 3 
bit phase shifters. Figures 9 through 12 illustrate the range 
of sidelobes, null depth and null location that results from 
quantization errors. The statistics from a large number of 
calculated patterns at 61 scan angles (0° to 60° in 1° steps) 
are presented in tables 2 through 7. 

The running sum and symmetric running sum _ roundoff 
algorithms had the minimum pointing error, in most cases a 
factor of 10 smaller than the other two algorithms. All four 
algorithms were comparable as far as sidelobe levels are 


concerned. 
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Figure 9. 50 Element Array Bayliss Difference Beam Radiation 
Pattern for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 10. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Weighted Random Roundoff 
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Figure 11. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Running Sum Roundoff 
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Figure 12. 50 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Symmetric Running Sum Roundoff 
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TABLE 2. COMPARISON OF 26 ELEMENT ARRAYS WITH 25 dB SIDELOBE 
LEVEL USING DIFFERENT ALGORITHMS 
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TABLE 3. COMPARISON OF 26 ELEMENT ARRAYS WITH 40 dB SIDELOBE 
LEVEL USING DIFFERENT ALGORITHM 
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TABLE 4. COMPARISON OF 50 ELEMENT ARRAYS WITH 25 dB SIDELOBE 
LEVEL ARRAY USING DIFFERENT ALGORITHMS 





TABLE 5. COMPARISON OF 50 ELEMENT ARRAYS WITH 40 dB SIDELOBE 
LEVEL ARRAY USING DIFFERENT ALGORITHM 
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TABLE 6. COMPARISON OF 76 ELEMENT ARRAYS WITH 25 dB DESIRED 
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TABLE 7. COMPARISON OF 76 ELEMENT ARRAYS WITH 40 dB SIDELOBE 
LEVEL USING DIFFERENT ALGORITHM 
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B. COMPARISON BY PHASE SHIFTER BITSIZE 

In this section, the phase shifter bit size is varied 
while all of the other array parameters are held fixed. For 
all algorithms increasing the number of bits improved the 
Sidelobe level, null depth, and beam pointing error. Therefore 
the data presented in this section will be limited to the 
symmetric running sum algorithm because the trends for the 
other roundoff algorithms are similar. The data is summarized 


in Table 8. Figures 13, 14, and 15 show the radiation pattern 


TABLE 8. COMPARISON OF RADIATION CHARACTERISTICS OF ANTENNA 
ARRAYS USING SYMMETRIC RUNNING SUM ALGORITHM 




















Mier i ln ener er Dar era I i) aC nr rn oy epee Tie ee ometeerel! 5 cll beraen’ inn 













en sbe se | se et eMh os fen 6 + 610.4 0 6.5 0.6 6.6 60s 46 0 OMe oe te Ong 
Filganare oc +'e o.suheit 0 + +-0 sil ore 17) I macgcers’ 





Pat Seet ot 





+r eretetetetele ToT + tet 














$70. 
Pai OMRON St EE ae NOt OT EDEN Oe CO Nt RE Oe he hd he Oe be eater S 
*. sete. 





re*e"e* 
ooo 5 ho 8o tes te 


| | cam) | Gail negree) 
-19.09 / 4.56 
4 [40 |- 4.00 / 2.37 _|-87.43 | -20.76 / 5.63 
-25.79 / 3.42 
-16.18 / 6.48 
-21.21 / 2.54 
-24.14 / 1.05 
40 | 3.37 / 0.111 | -84.18 | -22.32 / 6.79 
, 
| 
| 
: 
: 






ley SF Pe IED at ae ar Noe ata eK at he Sta at a ar a ae ark ea aCe ae et | 
Se RED SIE eee te eee ar a te) 
. ooo 2 © 6 @ GS BOOED GOO tt S.C. tte CN ei eR a at ee ae Se ae 










a0 


“27.19 / 5.12 
=16.98 / 6.24 
-22.64 / 2.18 
-24.61 / 0.57 
=23.66 / 4.39 
-28.66 / 3.83 


4a 


rae: 


“IoInN [TN |FN PN So Foor po lm [NR TN [NM FN 
NIN TNT nD THD |O |O | O OTN TO TO TO 





30 


for a 76 element array using 3, 4, and 5 bit phase shifters, 


respectfully. 
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Figure 13. 76 Element Array Bayliss Difference Beam for a 3 
Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 14. 76 Element Array Bayliss Difference Beam Radiation 
Pattern for a 4 Bit Phase Shifter using Symmetric Running Sum 
Roundoff 
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Figure 15. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 


33 


Figures 16, 17, and 18 compare the pointing error, null 


depth, and sidelobe level as a function of bitsize. 
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Figure 16. Pointing Error versus Scan Angle for a 76 Element 
Array with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 17. Null Depth versus Scan Angle for a 76 Element Array 


with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 
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Figure 18. Maximum Sidelobe Level versus Scan Angle for a 76 


Element Array with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 


36 


C. COMPARISON BY THE NUMBER OF ARRAY ELEMENTS 

In this section, the number of array elements is varied 
for each bitsize and roundoff algorithm. Figures 19, 20, and 
21 are plots of 26, 50, and 76 element arrays with a 3 bit 
phase shifter using the running sum roundoff algorithm. The 
results are similar to those in the previous section. It can 
be seen in these figures, that as the number of radiating 
elements increased, the sidelobe level increase above the 
desired level is smaller, as expected. A larger sampling of 
the various combinations of the number of radiating elements, 
bitsize of the phase shifter, and type of roundoff algorithm 
1s included in Appendix A. Appendix B has plots of null depth, 
pointing error, and sidelobe level as a function of scan 


angle. 


D. FORMULAS FOR EVALUATION OF THE AVERAGE SIDELOBE LEVEL 

It has been shown that the maximum sidelobe level of an 
antenna array is determined by the variance of the 
quantization error [Ref. 3: p. 2]. By taking the expected 
value of the radiation pattern, an expression for the mean 
error level is obtained. This was first done for reflectors by 
Ruze [Ref. 3]. A similar approach for discrete arrays of 


elements yields, 


D seme? LOS (8) , Pp (v-1) 
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Figure 19. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 20. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 21. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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wae re Ps is the expected value of the normalized power 


pattern (i.e., radiation pattern) and 


Wo is the variance of the quantization error 


6 is the pattern angle 
N 1s the number of radiating elements 
iS) is the error free power pattern 


The quantization error variance is calculated using, 


N 


oy Wa W) (V-2) 
= 
New 


where W is the quantization error at the nth element and wis 
the mean of the quantization error. 

Equation V-1 gives an estimate of the sidelobe level based 
on the variance of the quantization error resulting from a 
roundoff algorithm. Since the statistics of the quantization 
error for all the roundoff algorithms is determined primarily 
by the bitsize, the variances associated with each algorithm 
did not differ by large amounts. Equations V-1 and V-2 
explain why the sidelobe levels were essentially independent 
of roundoff algorithm. Figure 22 is a plot of the average 
sidelobe level predicted by the first term on the right hand 
Side of V-1 and a power pattern calculated using weighted 


random roundoff. 
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Figure 22. Plot of Predicted Average Sidelobe Level and Actual 
Power Pattern using Weighted Random Roundoff 
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E. FORMULAS FOR CALCULATING THE BEAM POINTING ERROR 

An accurate determination of target location by a radar is 
made possible by the difference beam null position. A closed 
form expression has been derived by Frank and Ruze [Ref. l:p. 
11-37,38] for the beam pointing error when regular roundoff is 


used. The difference pattern can be written as 





E(0) =23(sin(2%sin (6) -y,) +sin ( “FS sin (8) ~y.)+... 


(V-3) 
+sin (EES sin (8) - wn) ) 
2 


where Wf now includes the phase necessary to scan the beam to 
6 


weeds well as any quantization error. For a difference 


pattern at 0, the beam pointing error, 60 is 


80= sin (422°1) 85 sin (6) -y,) . (v-4) 
n=1 


If the phases are nearly those required to produce a null at 


0,, the sine can be replaced by its argument, and equation V-4 


reduces to, 


N 
2 
neseiher’ (Q) =a Se (V-5) 


2b t size-1 
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Figure 23 is a plot of the predicted pointing error using 
equation V-5 versus the calculated pointing errors for regular 


roundoff. 
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Figure 23. Predicted Pointing Error (*) and Calculated 
Pointing Error (+) versus Scan Angle for a 50 Element Array 
with a 4 Bit Phase Shifter using Regular Roundoff 
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V. CONCLUSIONS 

The running sum roundoff algorithms reduce the beam 
pointing error by a factor of 10 compared to the other 
roundoff algorithms, without degrading the maximum sidelobe 
levels or gain performance. The symmetric running sum roundoff 
algorithm does not seem to offer any significant performance 
advantage over the running sum algorithm. The only possible 
advantage the symmetric algorithm may offer is in beam 
scanning rate. Since the roundoff algorithm is only applied to 
N/2 elements as opposed to N elements, a faster beam scanning 
rate may result. 

Even though the running sum algorithm yielded essentially 
zero beam shift and 60 to 80 decibel null depths, these 
numbers will not be achieved in practice due to manufacturing 
and assembly errors. Phase shifters are designed to a set of 
electrical and mechanical tolerances. In a 4 bit phase 
shifter, for example, the least significant bit may not be 
exactly 22.5°. It could be slightly more or less depending on 
how that individual phase shifter’s circuits were fabricated 
and assembled. These errors tend to be random from phase 
shifter to phase shifter and therefore their effects can be 
estimated using statistical approaches such as those described 


in Chapter IV. 
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APPENDIX A —- PLOTS OF RADIATION PATTERNS AND QUANTIZATION 
ERRORS 


This section will allow the reader to examine the detailed 
effect of the number of elements, bitsize, and the roundoff 
algorithm on the radiation pattern and quantization error. The 
radiation patterns were arbitrarily chosen for a scan angle of 
5 degrees using 26, 50, and 76 elements with 3, 4, and 5 bit 
phase shifters using regular, weighted random, and running sum 


algorithms. The data is contained in figures 24 through 62. 
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Figure 24. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 25. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 26. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 27. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter uSing Running Sum Roundoff 
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Figure 28. 26 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetrical Running Sum 
Roundoff 
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Figure 29. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 30. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 31. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 
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Figure 32. 26 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 33. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 
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Figure 34. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 35. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 36. 26 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 37. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 38. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 39. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 40. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 41. 50 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetrical Running Sum 
Roundoff 
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Figure 42. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 43. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 44. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 
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Figure 45. 50 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 46. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 
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Figure 47. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 48. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 49. 50 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 50. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter with Perfect Phase 
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Figure 51. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Regular Roundoff 
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Figure 52. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 53. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Running Sum Roundoff 
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Figure 54. 76 Element Array Difference Beam Radiation Pattern 
for a 3 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 55. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Regular Roundoff 
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Figure 56. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 57. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Running Sum Roundoff 


c> 









QUANT. ERROR (DEGREES) 
an 
c> 


79 


MAGNITUDE ,dB 





DEGREES 





50 

W 

Y 

O 

uw 

Y oLAA 4 [\ MA Af\ AI\ al A/\ Mi iWAW Ae 

y aA | VV WV VV" IV ayy (/V 

rf 

A 

re 

) es 

GC 80 

MAX SLL = -24.28 

PT. ERROR = 0.002 ——— NULL AT -84.85 


Figure 58. 76 Element Array Difference Beam Radiation Pattern 
for a 4 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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Figure 59. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Regular Roundoff 
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Figure 60. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Weighted Random Roundoff 
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Figure 61. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Running Sum Roundoff 
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Figure 62. 76 Element Array Difference Beam Radiation Pattern 
for a 5 Bit Phase Shifter using Symmetric Running Sum Roundoff 
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APPENDIX B —- RADIATION PATTERN CHARACTERISTICS 
This section will allow the reader to examine the 
radiation pattern characteristics for scan angle from 0 to 60 
degrees. The plots of the radiation pattern characteristics 
are grouped by array size, each using the same roundoff 
algorithm. In each figure the performance of 3,4, and 5 bit 
phase shifters is compared, in that order from top to bottom. 


The data is presented in figures 63 through 98. 
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Figure 63. Pointing Error versus Scan Angle for 26 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Regular 
Roundoff 
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Figure 64. Null Depth versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 65. Maximum Sidelobe Level versus Scan Angle for 26 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Regular Roundoff 
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Figure 66. Point Error versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifter using Weighted Random 
Roundoff 


88 


NULL DEPTH (dB) 


NULL DEPTH (dB) 





SCAN ANGLE (DEGREES) 


Lv 

aa 

= 

QO. 

LJ MAA av) 
[| 
aa 
S) 

> 

Zz 

-100 
0 10 


SCAN ANGLE —— 
Figure 67. Null Depth versus Scan Angle for 26 Element Arrays 


with 3, 4, and 5 Bit Phase Shifters using Weighted Random 
Roundoff 
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Figure 68. Maximum Sidelobe Level versus Scan Angle for 26 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Weighted Random Roundoff 
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Figure 69. Pointing Error versus Scan Angle for 26 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 
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Figure 70. Null Depth versus Scan Angle for 26 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 71. Maximum Sidelobe Level versus Scan Angle for 26 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Running Sum Roundoff 
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Figure 72. Pointing Error versus Scan Angle for 26 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 73. Null Depth versus Scan Angle for 26 Element Arrays 


with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 
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Figure 74. Maximum Sidelobe Level versus Scan Angle for 26 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 
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Figure 75. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4 and 5 Bit Phase Shifters using Regular 
Roundoff 
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Figure 76. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 77. Maximum Sidelobe Level versus Scan Angle for 50 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Regular Roundoff 
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Figure 78. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Weighted 
Random Roundoff 
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Figure 79. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Weighted Random 
Roundoff 
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Figure 80. Maximum Sidelobe Level versus Scan Angle for 50 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Weighted Random Roundoff 
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Figure 81. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 
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Figure 82. Null Depth versus Scan Angle for 50 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 83. Maximum Sidelobe Level versus Scan Angle for 50 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Running Sum Roundoff 


POINTING ERROR (DEGREES) 





POINTING ERROR (DEGREES) 


SCAN ANGLE (DEGREES) 


0.1 


=a 
0 10 20 30 40 90 60 


SCAN ANGLE (DEGREES) 


Figure 84. Pointing Error versus Scan Angle for 50 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 85. Null Depth versus Scan Angle for 50 Element Arrays 


with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sun Roundoff 
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Figure 86. Maximum Sidelobe Level versus Scan Angle for 50 


Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 


108 


POINTING ERROR (DEGREES) 


POINTING ERROR (DEGREES) 





POINTING ERROR (DEGREES) 


SCAN ANGLE (DEGREES) 


Figure 87. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Regular 
Roundoff 


NULL DEPTH (dB) 


NULL DEPTH (dB) 


NULL DEPTH (dB) 





SCAN ANGLE (DEGREES) 


Figure 88. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Regular Roundoff 
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Figure 90. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Weighted 
Random Roundoff 
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Figure 91. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Weighted Random 
Roundoff 
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Figure 92. 
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Figure 93. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Running Sum 
Roundoff 
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Figure 94. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Running Sum Roundoff 
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Figure 96. Pointing Error versus Scan Angle for 76 Element 
Arrays with 3, 4, and 5 Bit Phase Shifters using Symmetric 
Running Sum Roundoff 
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Figure 97. Null Depth versus Scan Angle for 76 Element Arrays 
with 3, 4, and 5 Bit Phase Shifters using Symmetric Running 
Sum Roundoff 
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Figure 98. Maximum Sidelobe Level versus Scan Angle for 76 
Element Arrays with 3, 4, and 5 Bit Phase Shifters using 
Symmetric Running Sum Roundoff 
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APPENDIX C = COMPUTER PROGRAMS 
This section includes all the major programs used in 
compiling the data for this thesis. The programs for plotting 
and saving the data are not included in this appendix. All 


programs were written in MATLAB. 


A. MAIN PROGRAM 
Computes the Radiation pattern of a linear array with either 
a Bayliss, Taylor, or a Uniform Distribution using any size 


of Antenna Array with any Bitsize Phase Shifter. 


begin _scan_angle=0; Caintelalizing the anatial 
end _scan_angle=60; * conditions for the program 


Scart angle=-90; 
stop_angle=90; 

delta=.5; 

num_elements=[26 50]; 

phase algorithm=1; 

bit _size=(3 4 5]; 
selection=7; 
spacing per wavelength=.3; 
theta_step=1; 

Wl=1; w2=1; wW3=1; w4=1; w5=0; 


G1=0; d2=1;: da3=0; 


CALE 


rl=O0; r2=1; r3=0; r4=0; 

desired algorithms=[2]; 

sl1=25; 

nbar=5; 

elem counter=0; 

w6=[]; 

w7='no'; 

mainmenu % Calls a pulldown menu program for parameters 
% modifications 

elke 

bk=2*pi1; 

C=p1/ leo; 

del=.5; 


it=round((stop angle-start angle)/delta)+1; 


pe=0; 
while selection == 6 
al=length(num_elements) ; % The "a" varibles are used to 


a2=length(bit_ size) ; initialize the size of the 


ovo 


a3=length(desired_ algorithms) ;% arrays depending on the 


oe 


selections 

a4=al*a2*a3; % made for the particular 
% simulation 

a5=length(begin_scan_angle:theta_step:end_scan_angle) ; 


nel=max(num_elements) ; 


pterr=zeros(a4,a5); % Point Error Matrix 
Dred err=Zeres( 44 705)-- % Predicted Pointing Error Matrix 


eee 


mebldapth=zeros(a4,a5); $ Null Depth Matrix 


sllmax=zeros(a4,a5); % Maximum Sidelobe Level Matrix 
ang=zeros(a5,it); % Angle Matrix 

e=zeros(a5,it); % E Field Matrix 
phss=zeros(a5,16); % Phase State Matrix 
amp=zeros(a5,nel); %$ Element Amplitude Matrix 
XS1=zeros(a5,nel); % Phase Matrix 
qph=zeros(a5,nel); $ Quantized Phase Matrix 
trunc=zeros(a5,nel); % Quantization Error Matrix 
scan=zeros(1,a5); % Scan Angle Matrix 
angl=zeros(a5,2001); % Angle Matrix for -1 to +1 degrees 


% about Scan angle 
el=zeros(a5,2001); %$ E field Matrix for -1 to +1 
% degrees about Scan Angle 


count=0; 


fer eclement=num_ elements 
elem counter—-elem counter+1; 
for phasemod=desired algorithms 
for nbits=bit size 
count=counttl; 
ths_count=0; 
plot _counter=1; 
for ths=begin_scan_angle:theta_step:end_scan_angle 
eisecount—-ths counttl; 


minl=0; $ Initializing Variables to Determine 


23 


max l—=—120; $ Pointing Error, Null Depth, and 
sllmaxl=-120; % Maximum Sidelobe Level 
sllmax2=-120; 


Scan( ths COuUnEy) Seliar, 


% BAYLISS Coefficients for Difference Beams 


if dl ==1 $ Calling the Appropriate 
taylor % program for the Selected 
elseif d2 == 1 $ Distribution 
bayliss 
elseif d3 == 
uniform2 
end 


$¢ Determining the maximum Amplitude to Normalize the 
% Amplitudes 


am=max(abs(amp(ths count, 1l:num elements(elem counter) )jme 


if am == 
am =l; 
end 


amp(ths count, :)=amp(ths_ count,:)/am; 


% Scan Angles in Degrees 
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us=sin(ths*r); 


psis=bk*spacing per _wavelength*us; 


% Generate exact Phase Required at Each Element. When 

% PHASEMOD =0, the Phase is exact and the Phase Slope is 

$ Continuous ( no Modulo 360 or 2pi. Note that Positive 

% Scan Corresponds to Increasing Phase Lag with Increasing 


$n. 


for i=l:num_elements(elem_ counter) 

if sym == 1 & phasemod == 3 

xSi(ths_ count,1)=(-1)*(2*i-(num_elements(elem_counter)+1))/2 
*psis; 

else 

oui ths count, 1)=(-1)*(1i-1)*psis; 

end 


end 


ie XSi(ths count,1) 
< xsi(ths_count,num_elements (elem_counter) ) 
xmin=xsi(ths_count,1); 
else 
xmin=xsi(ths count,num_elements(elem_ counter) ); 


end 


%$ Normalize to the Most Negative Phase so that all xsi(i)> 0 
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xSi(ths count, :)=xsi(thiscoune, = ) jana, 


qph(ths_count,:)=xsi(ths_count,:); 


%* number of BITS is nbit; Number of Phase States is 


$ 2*(nbit); Each Bit is 360/(2*nbit) Degrees 


if phasemod > 0 
nstates=2“*(nbits); 


bitsize=2*pi/nstates; 


for i=l:nstates 


phss(ths_count,1)=bitsize*(i-1); 


end 
if phasemod == 1 ¢ Calling the Desired 
truncate * Roundoff Algorithm 
elseif phasemod == 2 
Heo) UL 8 (Clone 26 
elseif phasemod == 3 
runround 
elseif phasemod == 
other 
end 
end 


% Determining the Quantization Error 
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mame ( ths Count,:)—GQph(ths count, :)—xsi(ths_count,:); 


$ Begin Pattern Loop 
fiotdi—L: lt 
th=start angle+(i-1)*delta; 
ang(ths_count,1)=th; 
u=sin(th*r); 
psi=bk*spacing per wavelength*u; 
sum=0; 
$ sum to get Array Factor 
Zorenn=!t:num elements(elem counter ) 
Sepuase-—Ggph(ths count,nn); 
if (nn > (num_elements(elem_counter)/2)) & (d2 == 1) 
q phase=q _ phasetpi; 


end 


Ssum=sumtamp(ths_count,nn) *exp(j]*(psi*(nn—-1)+q_phase) ); 


end 


% Normalize the Pattern to the Number of Elements 
©y—-abs(sum) /num_elements(elem counter); 

e(ths_ count,1)=ev%2; 

e(ths count,i)=10*log10(e(ths count,1i)); 


end 


if d2 == 1 ¢ Refining the Radiation 


2-2 


nuldepth 
elseif dl ==1 | d3 ==1 
max lobe 


end 


Peer (COUNT, CHSeecount)| peer, 


% increment to .001 Degrees 
$ to determine the Point 
$ Error, Null Depth, and 


% Maximum Sidelobe Level 


% Saving the Pointing Error 


nulldpth(count,ths count)=minl-maxl; % and Null Depth to the 


err pred 


plot nol 


end 


plot_no0d 


mat save 


$ appropriate Matrix 


% Determining the Predicted 


$ Pointing Error 


Calling the Plotting Program 
for the Radiation Pattern and 


Quantization Error 


Calling the Plotting Program 
for the predicted versus 


calculated pointing Error 


Saving the data to plot the 
various combinations of 
pointing error, null depth, and 
Max. Sidelobe level vs. Scan 


Angle 
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mlat(0:60,pterr(count,:)), grid, xlabel(’SCAN ANGLE(DEG)’), 

ylabel( ‘POINTING ERROR (DEGREES) ‘),title(labell) 

meta scantot 

plot(0:60,nulldpth(count,:)),grid,xlabel( ‘SCAN ANGLE(DEG)’), 
ylabel( ‘NULL DEPTH (dB)’),title(labell) 

meta scantot 

plot(0:60,sllmax(count,:)),grid,xlabel(’SCAN ANGLE (DEG)’'), 

ylabel( ‘MAXIMUM SLL (dB)’),title(labell) 

meta scantot 

!gpp scantot 


{del scantot.met 

savefile % Saving the data 
end 

end 

end 

if w3 == 

plotcomp % Plotting the the data save by "matsave" 


end 


selection =7; 


end 


Mya) 


B. DETERMINATION OF THE POINTING ERROR 
This program determines the pointing error, null depth, and 
maximum sidelobe level for a Taylor Distribution and 


normalizes the radiation pattern to 0 decibels. 


del=.5; 

1€2=—-1* (sStarerangle+), 
1t3=-1*(start_angle-1); 

for 1=((it2+ths)/del):.002:((it3+ths)/del) 
th=start_angle+i*del; 

angl(ths count, (i=((it2tehs) /del)) 50078) air 
u=sin(th*r); 

psi=bk*spacing per wavelength*u; 


sum=0; 


%$ sum to get Array Factor 

for n=l: numeelements(ellensecounter, 
q phase=gph(ths count,n),; 
if (n > (num elements(elem icotnter)72)) 
q phase=q_ phasetpi; 


end 


sum=sumt+tamp(ths_ count,n)*exp(j*(psi*(n-1)+q_ phase) ); 


end 


% Normalize the Pattern to the Number of Elements 
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evl=abs(sum)/num_elements(elem_ counter); 


dummy=evl*evl; 


if abs(dummy) < le-9 
dummy=le-9; 


end 


Sweeems Count, (1-((1t2+ths) /del) )*500+1)=10*log10(dummy) ; 
evl=0; 


dummy=0; 


$ Determining the pointing error and null depth 


if (el(ths_ count, (1-((1it2+ths)/del))*500+1)) < (minl) 
minl=el(ths count, (1-((1t2+ths)/del) )*500+1); 
theta=th; 


pe=theta-ths; 


end 

end 

mom a—1l:it $ Determining the location 
if e(ths count,i) > maxl %¢ of the first maximum 
Maxi=-e( ths Count, i); %* of the differenece beam 
thetamax=ang(ths_ count,i); $ radiation pattern 
k=1} 
end 


sie al 


end 


m=k; 


max2=maxl; 


%$ Determining the location of max null if the maximum of the 


$ radiation pattern lies on the positive side of the null 


if thetamax > theta 
£fOR21—le ee 
if e(ths count,1i) <= max2 
max2=e(ths count,1); 
k1=i; 
else 
break 
end 


end 


% Determining the location of the second maximum of the 


% difference beam radiation pattern 


Forres) 
if abs(max2) >= abs(e(ths_count,1)) 
max2=e(ths count,1); 


k2=1; 


IL 8) 


else 
break 
end 
end 
b1=k2; 


end 


$ Determining the location of maximum null if the max of the 


$ radiation pattern lies on the negative side of the null 


if theta > thetamax 


for i=m:it 
if e(ths_ count,i) <= max2 
max2=e(ths_count,1); 
k1l=i; 
else 
break 
end 


end 


$ Determining the location of the second maximum of the 


$ difference beam radiation pattern 


for i=kl:it 


IE S38: 


if abs(max2) >= abs(e(ths_count,i) ) 


max2=e(ths count,i); 


k3=1; 
else 
break 
end 
end 

b2=k3; 

end 

amax=0; 


for 1=(b1-2)- 1 
Pe (thsecOunes) aecmeclinax 


amax=e(ths_count,i); 


135. 
else 
break 
end 
end 
amax=0; 


FOr 1=(D242 1 ee 

if e(ths_count,i) < amax 

amax=e(ths_ count,i); 
11=1; 


else 
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Determining the location 
of the first null on the 
negative side of the scan 


angle 


Determining the location 
of the first null on the 
positive side of the scan 


angle 


break 


end 

end 

fem i-1:1 % Determining the maximum 
mf e(ths count,i) > sllmaxl %* sidelobe level on the 
sllimaxl=e(ths_count,1); % negative side of the 
end $ scan angle 

end 

fem I=ll:it % Determining the maximum 
if e(ths_count,i) > sllmax2 %¢ sidelobe level; on the 
sllmax2=e(ths_ count,1); % positive side of the 
end %$ scan angle 

end 


if (abs(amax)+3) >= abs(sllmax2) 


sllmax2=-1000; 


end 

if sllmaxl > sllmax2 % Determining the maximum 
sllmax(count,ths_count)=sllmaxl-maxl; % sidelobe level 

2 LES 


sllmax(count,ths count)=sllmax2-max1; 


end 


5 


e(ths_ count, :)=e(ths count, :)-maxl; 4 Normalize ene 
$ Radiation 
yb=min(el(ths count,:)); $ pattern 
if yb <= -89.90 
[xa,xb])=min(el(ths count,:)); 
el(count,xb)=-90+max1; 
minl=el(ths_count,xb); 
end 


el(ths count, :)=el(ths count, :)—-maxl]; 


X=((1lE2Ztthsiydel ja. 


y=((it3+ths)/del)+1; 


C. BAYLISS DISTRIBUTION 

Computes Linear BAYLISS Coefficients for Difference Beams 
mu=zeros(1,11); 

z=zeros(1,10); 


b=zeros(1,10); 


$ Linear Bayliss Distribution 


nbarl=nbar-1; 


fOr vi—i- 


mu(i)=1-.5; 


end 
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=e 25=-sll-25; 
z(1)=1.87+S1125* .038; 
mee )—2.5+s1125* .016; 
Zeon —3.35+Ss1125*.019; 
z(4)=4.25+51125*.016; 


a=1.45+s1125* .042; 


for ns=5:nbar 


Z(ns)=sqrt(a*atns%*2); 


end 


Sigma=mu(nbar+l)/z(nbar); 


for mms=l:nbar 


bb=1; 
for ns=l:nbarl 
bb=bb* (1-(mu(mms)/(Sigma*z(ns)))*2); 
end 

bbb=1; 


for lls=l:nbar 
DEG@emlssmnc yea O ye (¢(ils—mms) > 0) 
bbb=bbb* (1-(mu(mms)/mu(lls))*2); 


end 


end 


So) 


bes=(-1)*(mms); 

b(mms )=mu(mms)*2/bes*bb/bbb; 

end 

for alf=1: (num elements (elem countess), 2) 


rho=(2*alf—-1)/num elements(elem counter); 


for lis=l1:nbar 

pmu=mu(lls)*pi*rho; 

bes=sin(pmu); 

gg=ggtb(lls)*bes; 

end 
amp(ths_count,num_elements(elem_ counter) /2+1l-—alt )=abs (gene 
amp(ths_count,num_elements(elem_ counter)/2talf)=abs(gqg) ; 
end 


end 


D. ROUNDOFF ALGORITHMS 

1. Regular Roundoff 
% Computes the Phase at each Array element by Truncation 
for i=l:num_elements(elem_ counter) 
nsteps=round(xsi(ths_count,i)/bitsize); 
xs=nsteps*bitsize; 
qph(ths_count,1)=xs; 


end 
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2. Weighted Regular Roundoff 

% Program Quantizes Phase Slope using Random Roundoff 

for i=l:num_elements(elem_counter) 

nsteps=fix(xsi(ths_count,1)/bitsize) ; 
plo=nsteps*bitsize; 

phi=plo+bitsize; 

gpn¢ths count, 1)=plo; 

xn=rand; 

%$ Criteria for Weighted Random Roundoff 
phnorm=(xsi(ths count,1)-plo)/bitsize; 
rnd=pnorm+xn; 

sae foe etel emily, 
gph(ths_count,1i)=phi; 
end 
end 
3. Running Sum Roundoff 
% Computes the Phase at each Array element by “the Running 


Sum Method" and "“Symetric Running Sum Method" 


nloop=num_elements(elem_counter)/2; 
erse 
nloop=num_elements(elem_counter); 


end 


3s 9 


fOr 1-1: noes 
nsteps=fix(xsi(ths_count,1i)/bitsize) ; 
phi _lo=nsteps*bitsize; 


phi hi=phi_ lot+bitsize; 


sone=st+phi lo-xsi(ths_count,i); 

stwo=stphi hi-xsi(ths_count,i); 

if abs(stwo) < abs(sone) 

S=stwo; 

gph (ths count;1)—pha has, 

if s = 
gph(ths_count,num_ elements (elem_counter)+1—-1)=—-phi_hi; 
end 

end 

1f abs(sone) < abs(stwo) 

s=sone; 

gph(ths count,i)=phi_ lo; 

if sym == 
gph(ths_count,num elements(elem counter)+1-1)=-phi_1lo; 
end 

end 


end 
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